Effects of shadowing in oblique-incidence metal (100) epitaxial growth 
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The effects of shadowing in oblique incidence metal (100) epitaxial growth are studied using a 
simplified model. We find that many of the features observed in Cu(100) growth, including the 
existence of a transition from anisotropic mounds to ripples perpendicular to the beam, can be 
explained purely by geometrical effects. We also show that the formation of (111) facets is crucial 
to the development of ripples at large angles of incidence. A second transition to 'rods' with (111) 
facets oriented parallel to the beam is also found at high deposition angles and film thicknesses. 
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While a variety of surface relaxation processes, such as 
adatom diffusion on terraces and near steps, as well as 
edge-diffusion and detachment, are usually assumed to 
determine the surface morphology in epitaxial growth, 
recently it has been showni2i2ii that the deposition pro- 
cess can also play an important role. For example, 
in the case of unstable metal epitaxial growth with an 
Ehrlich-Schwoebel (ES) barrier 5 to diffusion over de- 
scending steps, the short-range (SR) attraction of de- 
positing atoms to the substrate can significantly enhance 
the surface roughness and selected mound angle. 4 In ad- 
dition, for sufficiently large angles of incidence, the de- 
position angle can also play an important role. For ex- 
ample, in recent experiments on Cu/Cu(100) epitaxial 
growth 1,2 at 250 K a gradual transition was observed 
from symmetric mound structures for deposition angles 
up to = 55° (where is the angle between the beam 
and the substrate normal), to asymmetric mounds with 
increasing slopes for deposition angles up to 70°, to asym- 
metric ripples oriented perpendicular to the beam with 
(113) / (111) facets on the shadow/illuminated sides at 
= 80°. Similar results have been obtained in grazing in- 
cidence Co/Cu(001) growth^ for which it was found that 
the resulting surface anisotropy also leads to strong uni- 
axial magnetic anisotropy. Thus, understanding the ef- 
fects of oblique incidence deposition is important since it 
may lead to the possibility of controlling both the surface 
morphology and magnetic properties in epitaxial growth. 

Although the effects of shadowing on thin-film mor- 
phology have been extensively studied in the case of 
amorphous and polycrystalline columnar growth^ the 
case of epitaxial growth is not as well understood. We 
note that in Ref. U it was shown that at large angles of 
incidence the long-range (LR) van der Waals attraction 
of depositing atoms to the substrate may play an impor- 
tant role in addition to the SR attraction. Therefore, 
a fully realistic simulation of oblique-incidence epitaxial 
growth can be very time-consuming since it must take 
into account both the LR and SR interactions. As a 
result, recent simulations have focused on the submono- 
layer regime 8 or at most on the first few layers of vicinal 
growth. 9 However, at high angles of incidence the effects 
of shadowing and crystal geometry also play an increas- 



ingly important role. Therefore, it is of interest to deter- 
mine to what extent these purely geometric effects may 
determine the surface morphology. 

Here we present the results of simulations carried out 
using a simplified model of fcc(100) epitaxial growth in 
which the effects of shadowing are included but not the 
additional modifying effects of the SR and LR attrac- 
tion. Our results indicate that many of the qualita- 
tive and semi-quantitative features observed in Cu(100) 
growth^ 2 including the existence of anisotropy in the 
submonolayer regime, can be explained purely by geo- 
metrical (shadowing) effects. In addition, we find that 
the formation of (111) facets is crucial to the formation 
of ripple structures at large angles of incidence. The de- 
pendence of the critical thickness for ripple formation on 
deposition angle is also studied and good scaling behav- 
ior is obtained for the correlation length perpendicular 
to the beam as a function of deposition angle and film 
thickness. At higher thicknesses we also observe a second 
transition to 'rods' with (111) facets oriented parallel to 
the beam. We also find that the surface width increases 
exponentially with deposition angle and exhibits excel- 
lent scaling as a function of film thickness. 

Except for deposition, our model is very similar to pre- 
vious model a 10 ! 11 used to study metal (100) growth at 
normal incidence in which the correct crystal geometry 
has been taken into account. In particular, atoms are de- 
posited with a (per site) deposition rate F, while adatoms 
(monomers) on a flat terrace are assumed to diffuse with 
hopping rate D. Since the ES barrier typically plays an 
important role in metal epitaxial growth, the rate for 
an adatom at a descending step-edge to diffuse over the 
step is given by Des = De~ EES ^ kBT where Ees is the 
Ehrlich-Schwoebel barrier. Compact islands are also as- 
sumed and accordingly a moderate amount of edge- and 
corner-diffusion is also included, while the attachment of 
atoms to existing islands is assumed to be irreversible. 
Thus the most important parameters in our model are 
the deposition angle the ratio D/F of the monomer 
diffusion rate to the deposition rate, and the magnitude 
of the ES barrier. 

To take into account the effects of shadowing and crys- 
tal geometry, each deposited atom is assumed to travel 
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FIG. 1: (a) Surface roughness (in ML) as function of film 
thickness for — 0° — 85° for case of slow deposition and 
moderate ES barrier (b) Relative width deviation Sw/w(t,0) 
as function of scaled deposition angle Ot 1 for = 10° — 88° 
and t — 5 — 100 ML for three different cases: fast deposition, 
high ES barrier (top), fast deposition moderate ES barrier 
(middle), and slow deposition, moderate ES barrier (bottom). 



ballistically until its distance to the closest substrate 
atom is less than or equal to the nearest-neighbor dis- 
tance. The depositing atom then 'cascades' randomly 
via downward funneling (DF) 12 from a site correspond- 
ing to this atom until it reaches a four- fold hollow site. 
Thus, in our model atoms deposited on (111) facets are 
assumed to diffuse essentially instantaneously to the ter- 
race below. Given the extremely low barriers for diffu- 
sion on metal (111) surfaces (approximately 0.05 eV for 
Cu(lll)ii) for temperatures which are not too low this 
is a very reasonable approximation. 

In most of our simulations two different deposition 
rates were used - one corresponding to a 'slow' depo- 
sition rate (D/F = 10 5 ) and the other corresponding to 
a 'fast' deposition rate (D/F = 5000). Similarly, two 
different values of the ES barrier were also used, one cor- 
responding to a 'moderate' barrier at room temperature 
(Ees = 0.07 eV) and the other corresponding to a large, 
effectively infinite ES barrier. As in the experiments of 
van Dijken et aUi£ the overall deposition rate was as- 
sumed to be independent of deposition angle. Similarly, 
the azimuthal angle was chosen such that the deposition 
direction was parallel to the close-packed step-edge, i.e. 
along the [110] direction. 

In order to understand the dependence of the surface 



morphology on deposition conditions we have calculated 
a variety of different quantities as a function of average 
film thickness t (where t is in ML) and deposition angle 
0. These include the r.m.s. surface height or "width" w, 
the lateral correlation lengths £m and £j_ determined from 
the zero-crossing of the height-height correlation func- 
tions parallel and perpendicular to the beam, and the 
anisotropy a = Our simulations were typically 

carried out using relatively large system sizes (L = 512) 
and averaged over 100 runs. 

Fig. 1(a) shows typical results for the surface width as 
a function of film thickness for deposition angles ranging 
from 0° (normal incidence) to 85°, for the case of slow 
deposition and a moderate ES barrier. As in the Cu(100) 
growth experiments of van Dijken et al, 1 for < 50° the 
effects of oblique incidence and shadowing on the surface 
roughness are relatively weak. However, for larger depo- 
sition angles the surface roughness increases significantly 
with increasing deposition angle. As a result, the value 
of the effective roughening exponent j3 at large film thick- 
ness t (where w ~ t@) increases from (3 ~ 1/4 for small 
angles to a value close to 1 at = 85°. 

Fig. 1(b) shows the corresponding results for the rela- 
tive deviation in the width w(t, 9) compared to the width 
at normal incidence as a function of the scaling variable 
Ot 1 where the scaling exponent 7 ~ 0.075. Also in- 
cluded are similar results obtained for the cases of fast 
deposition with a moderate ES barrier, and fast depo- 
sition with a large ES barrier, which have been shifted 
up by factors of 10 and 100 respectively for clarity. As 
can be seen, there is excellent scaling for all three sets 
of deposition conditions. In addition, we note that the 
(unshifted) scaling functions in all three cases are almost 
identical except for some small deviations for large Ot 1 . 
The linearity of the scaling function in each case over 
approximately 4 decades shows that for fixed deposition 
rate the relative width deviation increases exponentially 
with increasing deposition angle 0. The small value of the 
scaling exponent 7 also indicates that the surface width 
depends much more strongly on deposition angle than on 
film thickness. 

Fig. 2 shows typical pictures of the surface morphol- 
ogy at a film thickness of 50 ML for different deposition 
angles for the case of slow deposition (moderate ES bar- 
rier). At a deposition angle ~ 70° (Fig. 2(a)) shad- 
owing leads to a preference for mounds to coalesce along 
the direction perpendicular to the beam and the onset of 
ripple formation. At a somewhat larger angle (0 c± 80°) 
asymmetric ripples with (111) facets on the illuminated 
side and (113) facets on the shadow side are formed (see 
Fig. 2(b) and inset) as in the experiments of Dijken et 
ali 2 The asymmetry of these ripples can be more clearly 
seen at c± 85° (Fig. 2(c)) which shows clearly the (111) 
facets on the illuminated side. We note that at this an- 
gle there is also evidence for a competition between the 
growth of asymmetric ripples perpendicular to the beam 
direction, and the growth of 'rods' with (111) facets par- 
allel to the beam. In simulations with reduced corner 
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FIG. 2: Gray-scale pictures of surface morphology (system 
size L — 512) for same conditions as in Fig. 1(a) at t — 50 
ML. Arrow indicates deposition direction. Inset shows ripple 
profile taken along line in (b). 



diffusion this competition leads to an extended regime 
in which isotropic pyramidal structures with (111) facets 
on all four sides are formed, similar to what is observed 
experimentally in Ref. 1. Finally, at = 88° (Fig. 2(d)) 
the perpendicular ripples are completely replaced by the 
growth of 'rods' with (111) facets growing parallel to the 
beam. 

The dependence of the anisotropy on both deposition 
angle and film thickness for this set of deposition con- 
ditions can be seen more quantitatively in Fig. 3. At 
a critical thickness t c , which increases rapidly with de- 
creasing deposition angle, there is an abrupt increase 
in the anisotropy, corresponding to the onset of ripple 
formation. The anisotropy then saturates as the rip- 
ples continue to elongate as well as coarsen but then de- 
creases as t _1//2 at higher film thicknesses. As can be 
seen in Fig. 3(b), at large angles there is already a small 
but noticeable anisotropy in the sub monolayer regime, 
even in the absence of steering effects due to SR and 
LR attraction™ This anisotropy is primarily due to the 
fact that shadowing tends to inhibit (enhance) the coa- 
lescence of islands along the directions parallel (perpen- 
dicular) to the beam. 

A detailed examination of the surface morphology in 
the multilayer regime indicates that the abrupt increase 
of the anisotropy and formation of ripples is due to the 
formation of (111) facets on the illuminated sides of 
mounds. Since (111) facets can efficiently capture and 
transport depositing atoms to the sides, this leads to a 
strong enhancement of mound coalescence in the direc- 
tion perpendicular to the beam, followed by the forma- 
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FIG. 3: Anisotropy a = as a function of film thickness 

for the case of slow deposition with a moderate ES barrier, 
(a) multilayer regime (b) submonolayer regime. In (b) the 
corresponding results for a large ES barrier are also shown 
(open symbols). 



tion and growth of ripples with extended (111) facets on 
the illuminated side. However, the saturation and even- 
tual decrease of the anisotropy may be attributed to the 
formation of structures with (111) facets on the sides par- 
allel to the beam as well as on the illuminated side (see 
Fig. 2(d)) whose growth competes with ripple growth. 
Such structures are particularly stable at large deposi- 
tion angle, while due to the large flux on the illuminated 
facet they tend to grow linearly in the direction of the 
beam. The decay of the anisotropy with exponent —1/2 
(see Fig. 3(a)) may then be explained by the fact that in 
this regime the correlation length perpendicular to the 
beam grows as ^_\_ ~ t 1 / 2 while the correlation length 
parallel to the beam grows linearly with film thickness, 
e.g. f || - t. 

Fig. 4 shows the scaled perpendicular correlation 
length as a function of scaled deposition angle Ot 21 where 
7 ~ 0.075, for t = 5 - 100 ML and = 55° - 85° for 
three different sets of deposition parameters (see cap- 
tion) where (a) and (b) have been shifted upwards by 
factors of 100 and 10 for clarity. As can be seen, in each 
case there is excellent scaling over a wide range of film 
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FIG. 4: Scaled correlation length £±/£ 27 as function of scaled 
deposition angle for three different sets of deposition parame- 
ters: (a) fast deposition, large ES barrier (b) fast deposition, 
moderate ES barrier, and (c) slow deposition, moderate ES 
barrier. Slope of dashed lines is approximately 2.6. 

thicknesses and deposition angles. This scaling behavior 
also indicates that ripple formation may be expected to 
occur for deposition angles as low as 55°, although the 
relatively low value of 7 indicates that the critical thick- 
ness for ripple formation increases rapidly with decreas- 
ing deposition angle. We note that the slope of the scal- 
ing function after the onset of ripple formation (dashed 
lines in Fig. 4) corresponds to the t 1 ! 2 growth of rip- 
ples with (111) facets perpendicular to the beam, and is 
consistent with domain growth with non-conserved order 
parameter, 15 while the later 'saturation' corresponds to 
the slower coarsening of structures with (111) facets par- 
allel to the beam, and is consistent with one-dimensional 
fluctuation induced coarsening^ We also find excellent 
scaling for the scaled correlation length £j_/7d as a func- 
tion of the scaled film thickness tjl\ where Id — (D/F) 1 ^ 6 
is the diffusion length, although the resulting scaling 
function depends on deposition angle. 

Finally, we discuss the factors that may affect ripple 
formation in oblique-incidence epitaxial growth. Since 
ripple formation is associated with the formation of (111) 
facets on the illuminated sides of mounds, we expect that 
anything that promotes this tends to enhance ripple for- 
mation. For example, since a large ES barrier leads to 



faster mound formation as well as a larger selected mound 
slope, increasing the ES barrier leads to earlier ripple 
formation as shown in Fig. 4. Similarly, decreasing the 
flux and thus increasing D/F, or increasing the rates 
of edge- and/or corner-diffusion delays the formation of 
(111) facets and ripples as again shown in Fig. 4. We 
have also considered the case of low-temperature growth 
(D/F ~0) but assuming that DF as well as fast diffusion 
on (111) facets remain active. Somewhat surprisingly, we 
find that even in this case, for which the surface is 'sta- 
ble' at normal incidence, 10 the formation of well-defined 
ripples with (111) facets is observed at large deposition 
angles (0 > 60°). Thus, we conclude that while ripple 
formation may be enhanced by the presence of a mound 
instability, it is primarily a geometric effect due to shad- 
owing as well as the existence of rapid diffusion on (111) 
facets. 

In conclusion, we have used a simplified model to study 
the evolution of the surface morphology and anisotropy 
in oblique-incidence metal (100) growth. Our results 
indicate that much of the qualitative and even semi- 
quantitative behavior can be explained by geometric ef- 
fects which dominate at large deposition angles. Our re- 
sults also indicate the existence of a competition between 
ripples perpendicular to the beam and structures with 
(111) facets parallel to the beam which eventually leads 
to a novel transition from perpendicular ripples to 'rods' 
at high deposition angles and film thicknesses. While we 
expect that the inclusion of SR and LR attraction will 
affect the results presented here, the qualitative picture is 
unlikely to change. In particular, due to the extreme sen- 
sitivity on deposition angle, the most significant effects of 
including such interactions are likely to be a slight shift 
in the angles at which ripple and/or rod formation take 
place as well as a slight modification of the anisotropy. 
It is also possible that the inclusion of the SR and LR 
attraction may further stabilize the "isotropic phase" cor- 
responding to the competition between rods and ripples 
as is observed experimentally. In future work, we plan 
to carry out multiscale simulations in order to determine 
in more detail to what extent the effects of SR and LR 
attraction may modify this picture. 
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